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Abstract 

We compute analytically the distribution function P(E) for the energy 
E acquired by a Fermi gas after being subjected to an arbitrary time- 
dependent external potential (switching event). We relate the distribution 
function to a solution of a matrix Riemann-Hilbert problem and present 
explicit formulae for the low order cumulants of P(E). These general 
results are used to find the distribution of dissipated energy in a biased 
quantum point contact. 

1 Introduction 

A rapid switching of an impurity potential in a Fermi gas causes creation of 
electron- hole pairs which carry away some energy (so called "shake-up" process). 
This non-trivial many body effect was first studied more than 30 years ago in 
relation to the A-ray absorption in metals H} [^]- In this case the perturbing 
potential appear instantly (i.e. the switching time T s is of the order of the inverse 
conduction band width D) and the distribution P{E) of energy E carried away 
by electron hole pairs is a power law P(E) ~ E~ a in the threshold region 
E < h/T s ~ D. 

It was later found that the same physics is important in the absorption of 
ions by metallic sur faces |, |. In this case the perturbing potential V(r, t) 
changes slowly (i.e. T s ^> H/D). The distribution P(E) was found theoretically 
for arbitrary energies |3|, || (the result is also quoted in Eq. (Q) below) provided 
that the potential V(r,t) only causes s-wave scattering. 

Recent advances in fabrication technology have made it possible to subject a 
Fermi gas to potentials with non-trivial time and space dependence. For example 
such potentials are used in adiabatic electron pumps (see e.g. @]). While the 
distribution of the pumped charge has received a lot of attention recently 0, 
the closely related distribution of dissipated energy is still unknown. In this case 
the time-dependent potential V(r, t) changes slowly in time (as in the case of ion 
absorption) but neither the results nor the methods from j| are applicable since 
other (non s-wave) scattering channels are important. Finding the distribution 
P{E) for an arbitrary (slow) time-dependent potential turns out to be a non- 
trivial problem and some novel ideas are required to solve it. 

We study the ideal Fermi gas in the presence of a time-dependent potential 
V = V(r, t) assuming that the potential is zero in the distant past (t — > — oo) and 
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becomes time independent in the distant future. We compute the probability 
P{E) that the gas acquires energy E after this switching process. The results 
are expressed through the scattering matrix S(Ep, t) = S(t) at the Fermi energy 
Ef on the instantaneous value of the potential V(r,i). It is more convenient to 
deal with the characteristic function 

X (A) = J e- lEX P{E)dE (1) 

rather than with the probability distribution itself. In the rest of the paper 
we focus on the case when the scattering matrix changes slowly on the scale of 
Wigner's delay time, i.e. 

Here and everywhere below we use units with K = 1. We start by briefly 
summarising the facts which are already established in the literature (Eqs. ^ - 
0) below). 

At low energies E <C 1/T S the distribution function has the power law asymp- 
tote g§ 

P(E)=C(E-AEr, a= tr(ln^(co)) _ li (jJ) 

where AE is the difference between the ground state energies of the Fermi gas 
with and without potential V^r, oo) |]l0| and is given by 

i f° 

AE = — dE tr hi S(E,oo). (4) 

The average absorbed energy (E) = J P(E)EdE and the probability Poo to 
remain in the ground state are given by 

(E) = ~f dUx{(d t S)S^d t S)S^}, (5) 



Poo - expj^^ ix(D(u)D(-u))udu^ 



(6) 



respectively @, where D(u) = J e luJt In S(t)dt. When V(r,oo) ^ the inte- 
gral in (g) diverges logarithmically giving rise to the Anderson orthogonality 
catastrophe Jl^] . 

Finally, in the "commutative case" when the scattering matrix S(t) can be 
diagonalised in a time-independent basis (e.g. when there is only s-wave scat- 
tering) the problem admits a complete solution and the characteristic function 
is given by [§ |) 

ln X (A) = -iAEX + — / \D{uj)\ 2 u{1 - e- iuX )du. (7) 



4^ 2 Jo 

To illustrate the above results and compare them with our findings we con- 
sider a simple model of a quantum point contact biased by a rectangular voltage 
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pulse Vb(t) — Vo6(t)0(T s — t). After a gauge transformation (see Sec. ^ for de- 
tails) the potential becomes local in space and the scattering matrix takes the 
form 



S{t) 



Ae^W —B* 
B A*e- 1 ^ 



(8) 



where A, B are transmission and reflection amplitudes and 4>(t) — J_ eVf,(r)dr 
is the Faraday flux. 

The results (§-|) are of not much use for the simple model introduced above: 
since the matrices in (||) do not commute with each other at different times, 
Eq. (0) is not applicable and the only known result is Eq. (|J) for the average 
energy. This result, however, is just a straightforward combination of the Joule 
law (E) = J I(t)V b {t)dt with the Landauer-Biittiker formula @ [lj] for the 

current across the junction / = ^\A\ 2 Vb- 

In this paper we develop a technique for finding distribution function P{E) 
which allows us 

- to obtain the cumulants expansion of P(E) explicitly and to compute the sec- 



ond and third cumulants for an arbitrary switching process (see Eqs. (|2J, 25)) 

- to relate P(E) to a solution of a matrix Riemann-Hilbert (RH) problem 
(Eq. ©) 

- to find the distribution function of energy dissipation in the quantum point 
contact biased by a long rectangular pulse (see below) 

We finish the introduction by illustrating the relation between P(E) and the 
statistics of transmitted charge. Consider a quantum point contact described 
above. When the rectangular pulse is long i.e. T s ^> l/(eVo) we solve the RH 
problem asymptotically to logarithmic accuracy in n = cVqT s / '(2ir) ^ 1; this 
gives 

ln X (A) = nln{\A\ 2 e- leV « x + \B\ 2 ) + 

+ ^H^r)^(\A\ 2 e-^ + \B\ 2 ), (9) 

where we assume eV^X <C n. Since each electron reaching the contact contributes 
eV to the total dissipated energy if it is transmitted (probability \A\ 2 ) and 
contributes nothing if it is reflected (probability |B| 2 ) we expect at least for long 
times T s that characteristic function (^) can be obtained from the distribution 
of transmitted charge. The latter is known to be binomial |ll| 

B(N,k) = k] ^[ k)] \ A \ 2k \B\ 2{N - k) (10) 

where B(N, k) is the probability that out of N incident electrons k are trans- 
mitted. The distribution p(N) of the number of incident electrons can be found 
from the energy dissipation of a completely open contact (i.e. one with A = 1). 
The mean number of electrons transmitted through an open contact is given by 
(N) = n, and the distribution p(N) is Gaussian in the vicinity of its maximum 



(see Sec. 4.1 for details) 

p(N) - e -^(N-nf/{A\nn) &t n _ N „ ^ ( U ) 
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There are therefore two sources of shot noise in a quantum point contact - one 
is the fluctuations in the number of incident electrons and the other is the fluc- 
tuations in the number of reflected electrons. We are now able to establish that 
these two sources are statistically independent by observing that the character- 
istic function (|^) corresponds (with logarithmic accuracy) to the convolution of 
p(N) and B(N, k): 

P(E)=J2p(N)B(N,k)S{E~eV k). (12) 

N,k 

The physical picture of the two statistically independent sources of quantum 
shot noise was suggested (but not proven) in |l6|]. Note, that in order to confirm 
the statistical independence we need to obtain the sub-leading term in the large 
n expansion of x(A) (namely the term proportional to In n in Eq (^)) and analyse 
its dependence on the transmission amplitude A. We are not aware of any other 
technique capable of obtaining this term for a partially open channel. 

In fact, measuring the energy dissipated in a contact in the presence of a rect- 
angular bias pulse can be viewed as a "quantum charge measurement scheme" 
alternative to (and somewhat conceptually simpler than) using a precessing 
spin as was suggested in |l5| . The detailed comparison of the two measurement 
schemes will be the subject of future work. 



2 Overview of the method 

Let H(t) = Hq + V(r,t) be the time dependent Hamiltonian of the Fermi gas 
in the external potential, where Hq is the Hamiltonian of the free Fermi gas. 

If V(r, oo) 7^ then there are two contributions to the energy of the gas at 
t — > oo. The first is the sum of the energies of the electron-hole pairs and the 
second originates from the adiabatic shift of the single particle energy levels. If 
the energy levels in the vicinity of the Fermi surface shift uniformly by V(r, oo) 
then we can adiabatically switch off the potential V(r, oo) preserving the energy 
of existing electron-hole pairs and creating no additional ones. Therefore to 
account for V(r, oo) we should add —iAEX to the characteristic function x(A). 
We assume V(r, t > T s ) = in the rest of the paper. 

The evolution operator XJ t which relates the many-body wave function |\&(t)) 
to the the ground state wave function |\& ) via \^{t)) = Ut\^o) obeys the 
Schrodinger equation 

i^-Ut = H(t)U t . (13) 
at 

The characteristic function (|l]) is given by the average 

X (A) - lim (VolU^e-^Ute^Va). (14) 

t — >oc 

Usually such averages are evaluated by expanding both U and U~ l in the powers 
of V(r, t) and by evaluating the result using the Keldysh diagram technique (see 
for application of this method in the commutative case) . We follow a different 
approach suggested in || and separate the single-particle scattering problem 
from the averaging over the Fermi-gas ground state. The scattering problem 
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consists of finding the electron annihilation operator Ck(i) in the Heisenberg 
representation 

Pk(t) = Uf^Ut = e-^'^^^Cks (15) 

k' 

where £k is the energy of the electron with momentum k. The matrix <7kk' (t) 
becomes time independent at t > T s when the potential is switched off. Only 
this limiting value is needed for computation of the characteristic function ( |l4| ) . 
It can be expressed as a Fourier component of the scattering matrix 

a kk <M = J!T' I dtS nn ,(E,ty^-^ (16) 

r(f )v{E) j 

where v{E) is the density of states, d is the number of space dimensions, n = 
k/|k| and E — (ek + ek')/2- Eq. ( |l6| ) is the first term in the expansion in 
(S t f§S ,t §|) and is valid as long as (ek— ek'^ff < E Tne detailed calculation 
is presented in appendix [a]. 

Using ( |l5| ) , ([l6|) we can compute the limiting value of the Hamiltonian in the 
Heisenberg representation 

i?(too) = U t -jH U tx = Y, £ kc k c k + c l W ^c k/ , (17) 

k k,k' 

where 

HW ^fiS)/ iS,(B - i) T e ' ,,k ' , '' , '' i '- < 18 ' 

Now the characteristic function ( |l4| ) acquires the form of a time independent 
expectation value and can be evaluated using a linked-cluster expansion (see 
appendix [b| for details) 

ln X (A) = ^ J dtti{(d t S)S\d t S)&} + W (19) 

where W is the sum of the ring diagrams starting with the second order term 



2 W 3- 
Here the Green's function is 

G(p, At) = -ie ipAT (9(- P )9(At) - 9(p)6{-At)) (21) 
and the vertex is given by A(p — p', r) = M(p — p'){6{r) — 9(t — A) where 



M(p-p')=i J S\t)^e- l{p - p,)t dt. (22) 

We integrate over r and sum over matrix indices of M(p — p') in each vertex 
and sum over p in each line. 
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The diagram expansion ( |20| ) makes the computation of the ?i-th order cu- 
mulant 

\{ E ))=* ^ U=o (23) 

relatively straightforward. Indeed, since the n-th diagram is proportional to A™ 
at small A no more than n terms is needed for the n-th cumulant. This gives 
the following results for the second and the third cumulants 

1 f°° 

«^ 2 )) = ^2 J tv(M(p)M(- P ))pdp (24) 
1 [°° 

((^ 3 )) = 7^ min(p 1 ,p 2 )tr{M(p 2 -p 1 )[M(-p 2 ),M(p 1 )])dp 1 dp 2 + 



(2tt) 3 

1 
4^2 



ii{M{p)M{-p))p z dp. (25) 



Finding the characteristic function xW is a much more difficult problem 
which we relate to the RH problem described below. We start by introducing 
the matrix valued function S+(t,y) which arbitrary interpolates between 1 at 
y = and S(t + A) at y = 1 

S+(t,l)=S(t + X), S+(t,0) = l (26) 

and analogously the matrix valued function S—(t,y) interpolating between 1 
and S(t) 

S-{t,l) = S{t), SL(i,0) = l. (27) 

It is also convenient to define the new variables z — r + it and z = r — it. For a 
fixed < y < 1 the RH problem consists of finding two matrix valued functions 
f±(z,y) that satisfy the following conditions: 
(i) f+(z, y) is antianalytic when Re z > 
y) is antianalytic when Re z < 

(iii) /±(oo,y) = 1 

(iv) for real t 

/+(-** + 0, y)f-(-it - 0, yy 1 = Q(t, y), (28) 

where Q = S^S^. 

Once the solution of the RH problem (or approximation to it) is found 
the characteristic function of dissipated energy is given by the two-dimensional 
integral 

lnx(A) = ±- J ^{dS+S- 1 AdS-SZ 1 } 

- ihj triQ-'dQAdf^fl 1 } , (29) 

where A is the wedge product, the integrals are over the strip (0<y<l,— oo< 
t < +oo) and the orientation is chosen in such a way that J dy A dt is positive. 

The difficult step in the outlined procedure is the solution of the RH problem. 
We consider two special cases. In the quantum point contact case with a partic- 
ular ^-dependence of S'-matrix the RH problem can be solved asymptotically(see 
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section |4.2| ) . The "commutative" case when the S matrix can be diagonalised 
in some time-independent basis is discussed below. For pure s-wave scattering 
the S matrix has only one nontrivial eigenvalue S(t) = e ~ 2lS v(t) ^ where <5y(i) is 
the phase shift at the Fermi energy. The addition of more scattering channels 
is straightforward. Formulating (^8[ ^9j), we choose the interpolating functions 
S± in the form 



which gives Q 
the solution 



S-(t,y) = e - 2iySv ^\ S+ = e - 2 ^Mt+A) (30) 
-2iy(5 v {t)~8 v {t+\)) _ Employing the method from Q we obtain 



In /± 0,2/) 



n 



dt 



M«i) - 5 v {h + A) 
z + iti ± 



(31) 



Substituting @|o|) into ([|| gives (@). 

Equations (g§)7 @, © and Eq. 
results of the paper. 



from the introduction are the main 



3 Wess-Zumino action and the Riemann-Hilbert 
problem 

In this section we employ the Wess-Zumino action to find W given by the 
diagram series (^0|). Firstly, we show that the functional derivative of VV with 
respect to the vertex A 

SW 

j = 2i-— = -2iGAG - 2iGAGAG (32) 

SA y ' 

satisfies the equation for the current in the Wess-Zumino model. Observe that 
the Green's function G(p, t) given by Eq. ( ^l|) obeys the same equation as the 
Green's function of a massless left-moving particle in one dimension, p being 
the momentum of the particle and x is the conjugate variable to p 

(<9 r - d x )G(x -x',t- t') = -iS{x - x')5(t - t'). (33) 

Note, that in the coordinate representation the vertex A takes the form 

A(r ) x) = iS^(d(r)-e(r-X)). (34) 
ax 

Using ([}3|) we get for the n-th order [n > 2) of the perturbation expan- 
sion d§) 

(d T - d x )j n = -i[A,j n -i]. (35) 
The second order diagram can be presented in the form 

-\ J ^IL(p,u)A(p,u)A(-p,-w) 

where II(p, ui) — - - - ~'p S ^g n ^ is a polarisation operator. Varying this equa- 
tion over A we have the following equation for j'2 

(3r - d x )j 2 - —d x A. (36) 

7T 



7 



Combining (|35|) and fl36| ) we arrive at the implicit equation for W 

(d T -d x )j + i[A,j] = —d x A; J = 2<5¥. (37) 

This equation is very similar to the equation for currents in the Wess-Zumino 
model. To exploit this analogy we introduce the functional 

VV = W + — f tr( dgSS-^dx = W~ — I tr A 2 (x,T)dxdT (38) 
8tt J 8tt J 

which obeys the implicit equation 

-1 - SW 

(d T -d x )j + i[A,j] = —(d r + d x )A; j = 2i—. (39) 
2ir oA 

Polyakov and Wiegman(PW) solved ( |39| ) for the case when A has the 
form 

A = i(d T - d x )gg~ 1 \ g\ T ^oo, x ^oo = 1, (40) 
where g takes values in SU (N) . Under these conditions Eq. ([$9|) has the solution 

W(A) = -±J dxdTtrid^gg^d^gg- 1 } 

- ^ J dxdrdye^tYid^gg^dvgg-'dxgg- 1 }. (41) 

In the second term in ( |4l| ) the integration is over a membrane in the target 
space with boundary g(r, x). The current is given by j = j^(d T + d x )gg~ 1 . 

In our case the PW procedure is not applicable because the vertex A does 
not admit representation ([l0]). This can be demonstrated by using Eq. ( |32| ) 
for the current and observing that it does not decay at x, r — > oo, in con- 



tradiction to the PW prediction. For example, the second term in ( |32[) gives 



7a = 2i f i r; 7 7TT9 A(t' . x')dx'dr' which is not zero if x — > oo, 

J ^ J \x—x'+T—r' — ivsgn[r — r')\ z \ ' / ' 

t — > oo along the line x = —r. 

Boundary conditions that are more general than (|40| ) were considered in 
[[l9[ using a rather involved method. Instead, we overcome the problem 
by rotating the contour of integration over the energy e in diagram expan- 
sion (|20|) (the procedure is clearer in the energy representation) and take e = 
it; where £ is real and varies from — oo to oo. As the Fourier transform of 

[0(t) — 9(t — A)] is e j( e _ e /r 1 ; we also need to rotate A to keep the integrals 
convergent. In this way we obtain x(^) at A = — iot; where a is real and pos- 
itive. The characteristic function for a real A is obtained using analyticity of 
x(A) in the lower half-plane. The vertex in the coordinate representation [0 
is A(t, x) = iS^(x) 9S q^ {0(t + a) — 9(t)) and the rotated Green's function is 
G(x — x',t — t') = 2ir(T-T'-i{x-x')) 1 This Green's function satisfies equation 
2d z G = S(x — x')S(t — t') where we denote z — r + ix. 

Similar arguments as above lead to the equation for the current j — 2i^jj- 
(note that the current now goes to zero as z — > oo) 

2d z ~j-[A,j} = — d- z A. (42) 

7T 



Equation ([fOj) now takes the form 

2d z gg- 1 = A = iS \x)^±(9{t + a) - 6{r)) (43) 

and has the solution with g £ GL(N,C) satisfying g\ z ~toc = 1, 

9(r,x) = ( T< - Q ' T >° (44) 

Jy ' ; \ Sn(x)f(z), -a < t < 0. y J 

The conditions on the function /(f) = /(r — ix) are as follows 

(i) / is antianalytical in the complement of the two vertical (i.e. parallel to x 
axis) branch cuts at r = and t = —a 

(ii) /(oo) = 1 

(iii) Values of / on the opposite sides of the cuts are related by 

/(-0-a-ix) = S ,t (x)/(+0-a-ix) 
^(^/(-O-ia;) = /(+0-ia;). (45) 

Conditions (i)-(iii) describe a matrix Riemann-Hilbert problem. Once the de- 
composition ( ^3[) is found the solution of (jl^) is given by the Wess-Zumino 
functional 

= dxdrtrj^^- 1 ^^- 1 } 

- ^ J dxdTdye^tYid^gg^d^gg^dxgg- 1 }. (46) 

In the second integral g is continued on the membrane that has g(r,x) as a 
boundary. We choose the continuation g(r,x, y) with < y < 1 based on an 
arbitrary interpolation S(x,y) between S(x) = S(x, 1) and 1 = S(x,0) 

5(r '"' y) -\ SKx,y)f(z,y), -a<r<0. (47) 

The function f(z,y) has the same branch cuts and behavior at infinity as /(f) 
and obeys ( f45| ) at each value of y 

f(-0 - a- ix,y) = S~ 1 (x,y)f(+0 - a - ix,y) 
S-\x,y)f(-0-ix,y) = f(+0-ix,y). (48) 

Substituting (|47| ) into ( (46| ) and taking into account ([38|) we get (see appendix |^ 
for details) 

W = ~ [ tiidSS^dfi-a + Q-ix^jjf- 1 }-^- f tridSS^dfi-O-ix^)/- 1 }- 
- I dxtr{d x S(x)S- 1 (x)d x S(x)S^ 1 (x)}. (49) 

47T J 

where /3i is the plane r — —a and is the plane t = 0. Note that the last term 
in (|4^) exactly cancels the contribution from the first order diagram in (|l9|). 

The only remaining step is the analytic continuation of ( ^9| ) from real positive 
a to a = iX. The method of the continuation is illustrated on Fig. |. We let 
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X 



-a 



Figure 1: Branch cuts for the continuation of f(z,y) to the complex values of 



take the complex value a = a' + ia" and obtain the continuation of ( |47| ) in the 
form 



f(z,y), r<-a',T>0 
S- 1 (x + ^-r,y)f(z,y), -a'<r<0 



where 



(50) 



f(-0-a' -ix,y) = S~ l {x - a" ,y)f{+Q - a! - ix) 
S-\x,y)f(-0-ix,y) = f(+0-ix,y). (51) 

In the limiting case a = iX the surfaces r = —a' , r = coincide and ( [i"9|) 
takes the form (E9). 



4 Energy dissipation in a quantum point contact 



In this section we apply Eq. ( |29| ) to a biased quantum point contact. The contact 
is described by the Hamiltonian 



1 d 2 



(52) 



where the potential U(x) is local in space and therefore the scattering matrix 



A -B* 
B A* 



(53) 



is almost energy independent. The presence of the bias potential V& (t) is taken 
into account by adding the term eVb(t)0(x) to the Hamiltonian ([52]). The gauge 
transformation ip(x,t) — * ij:( y x,t)e^ 1 '^ 1 ^ 6 ^ where <p(t) — J*^ eV{T)dr makes 
the Hamiltonian local in space 



H(t) = U-i-^-5(x)<t>(t)) 



(54) 



with the scattering matrix given by (^ 
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In what follows, we consider the rectangular bias pulse with amplitude Vq 
and duration T s . The flux <p(t) is given by 



0, t < 

4>{t) = { eV t, 0<t<T s 
eV T s , t>T s . 



(55) 



The RH problem for the scattering matrix (J8|) can be solved in two limiting 
cases: when the contact is transparent (i.e. does not reflect electrons) and when 
the time T s is large. 



4.1 Transparent contact 

In the case of transparent contact the scattering matrix (0) acquires the form 



S(t) 



J<t>{t) o 



e -i<t>(t) 



(56) 



Since it remains diagonal at all times Eq. (]?]) can be used, which gives 



= -if 



2 r 



(A - T,)' In 



+ i -\ ('^') 2A 2 in^A 



(57) 



Classically all electrons in the energy strip (Ep,Ep — eVo) that can reach 
the contact are transmitted. The average number of such electrons is n — 
veV^vpTs = cVqTs I (2ir) , where vf is a Fermi velocity and v — \/{2ttvf) is the 
density of states per unit length. Since each transmitted electron contributes 
energy eVo to the total dissipated energy it is convenient to introduce a new 
rescaled distribution function 



P c (k) = eV P(eV k) (58) 

that measures energy in the units of eVb- The new distribution P c (k) depends 
only on the parameter n. We plot P c (k) for three different values of n in Fig |[ 
For small n the distribution function has a divergence at small k 

P c = Cn/k 1 - 2 " 2 at k < 1/n, (59) 

which is the Nozieres and De Dominicis asymptote (||) with the prefactor C n 
given by 

C n = 2nsin(27m 2 )e- 3 " 2 r(l - 2n 2 )(2 7 m) 2 " 2 - 1 . (60) 

In the other limit n ^> 1 the distribution looks like Gaussian with mean value n 
and the standard deviation a = ^\/lnn. The Gaussian approximation is valid 
only in the region k — n ~ 1 . The distribution function at k 3> n is determined 
by the term A 2 In A in (^) , which leads to the logarithmically divergent second 
cumulant. 
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Figure 2: P c (k) is the probability that the one-dimensional Fermi gas in the wire 
acquires the energy E — keVo after applying the voltage Vq during the time T s . 
The distribution P c (k) is parametrised by a single parameter n = eVoT s / '(2n), 
which is the expected number of electrons passed through the contact. 

(a) The Faraday flux <j>(t) as a function of time. The potential switches on at 
t = and switches off at t = T s . 

(b) P c (k) for n = 0.4 shows the Fermi-Edge singularity (the broken curve is the 
power law asymptote P c — C n /k 1 ~ 2n with C n given by Eq. (|6C|). 

(c) For n — 0.8 > l/v2 the behaviour at small k is also determined by the 
power law asymptote (broken curve). However, now the charge distribution is 
not divergent at small k. 

(d) , (e) P c (k) for the larger values of n crosses over to the Gaussian distribution 
(shown by the broken curve) with the mean value n and the standard deviation 
a = —\/\nn. 



4.2 Large T s asymptote of %(A) 



For a semi-transparent contact we can still solve Eq. (|28h asymptotically when 
T s > l/(eVo)- The matrix Q = S^(x + A, y)S(x, y) has the form 



Q = 



| A |2 e - 4 (0 + -0_) + | B | 2 -A*B*(e 
AB(e 



i<t>+ 



\A\ 2 e l( 



\B\ 2 



(61) 



where </>+ = y<p{x + A), 4>- = y4>{x). In the limit A <C T s it has the simple time 
dependence 

Q = Qo + Q+e* yeVox + Q_ e - iyeVox (62) 

where Qo is diagonal and Q-{Q+) are upper (lower) triangular matrices, which 
are almost time-independent. This suggest the following anzats for f± 

= h+ + f 1+ e- yeV(rz 



f- 



fo- + fi-e 



yeV z 



(63) 
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After substituting (|63|) into (pa) we get the approximate solution 



/+ = (l + Q + Q- 1 e-^)/ 0+ (64) 



/- = (1 - Qo X Q-e yeVo *)f Q - 
where /o+,/o- are diagonal matrices satisfying the equation 



with 



— Qo — Q+Qq 1( 5- — 



/o+/o_ — Qc 



|A| a e-*W+-*-) + |B| 2 





l 

| j4 |2 e --(* + -*-) + | B | 



The solution for ln/o_ is given by 

i , In f |^|2 e - l (0+(^)-0-(^')) + 

ln/o- = ^- / — * ; — Ux' 

2m J x — x + iv 



1 

-1 



(65) 



(66) 



(67) 



(68) 



Substituting /_ = A_/ -, where A- = 1 — Q 1 Q-, into j29| ) we get 

W = J ^{dS+S- 1 AdS^SZ 1 } - ^ J^ti{Q- 1 dQ hdA^AZ 1 } - 

- -p [ tr {Q^dQ A A-idlnfo^AZ 1 } . (69) 
4tt J d 



The first two integrals are proportional to T s , the last integral is proportional 
to ln(T s /A). After long but straightforward algebra we finally get (||). 



5 Conclusions 

We considered the energy dissipation in a Fermi gas subjected to a perturbation 
with time dependent scattering matrix of a general form. Result (|2^) relates 
the dissipated energy distribution function P(E) to a solution of the matrix 
Riemann-Hilbert problem (p8|). 

For a quantum point contact biased by a rectangular voltage pulse we com- 
puted P(E) analytically for long pulse durations and analysed the relation be- 
tween P(E) and distribution of transmitted charge. We also obtained second 
and third cumulants of the dissipated energy for a general switching process. 
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A One-particle scattering problem 

In this section we relate the electron annihilation operator c k in the Heisen- 
berg representation to the one-particle scattering matrix S(E, t) and derive 
Eqs. ( [l5[ |l6| ) from the overview. 

Solving the equation for the time evolution of c k 

Ck = -«£kCk - i ^2 V kk > (t)Ok' (70) 

k' 

perturbatively we get 

c k (t) = e-^lc^-if dt'Vwit'y^^'cv+i 2 [ dt>Vw(t>y^-^ 

\ J — oo J — oo 

rt' \ 

'^*( e k'- e k")t", 



dt"W(« / )e^'-^"^ c k " - ...J . 

In the limit t — > oo we get 

Ck(ioc) = e~ lek *- ^(<W - «F kk <(e k ,ek<))ck<- (72) 

k' 

Here we introduced the function -Fkk'( e > e') satisfying the integral equation 

/d 
^V kkl (e-e 1 )G kl (e 1 )F klk ,(ei,e'), (73) 

Kl ^ 

where V kk <(u)) = / <^Hck' (t)e _8w * and G k (e) = e _ e ^ + . Q is the retarded Green's 
function. To take into account the slowness of V kk /(t) we introduce the new 
variables E = (e + e')/2 and cj = e — e' and notice that F varies slowly 
when E changes but varies fast when u> changes. Using Wigner's representation 
Fkv(E,t) = J e~ iu ' t F kV (E,w)^ we get the approximate solution of © 



F kk .{E,t) = V kk/ (t)+J2v kkl G kl {E)F^(E,t). (74) 

ki 

The first order correction to is given by 

idVdG^ i^dGdF idVBF 
SF = 2~dt~dE + 2 V dEm - 2lK G dE- (?5) 
Eq. ( 74 ) is similar to the equation that one gets in quantum scattering the- 



ory (see [22 



jl30). 



Substituting F kk < given by ( |74| ) into (72) and neglecting the terms of the 
order §f f§ we get (|l6|). 

B Linked cluster expansion for xW 

From Eq. (|l^) we see that the characteristic function (jlj) is given by the linked 
cluster expansion 

ln X = ln(* |e-^ (oo) V^ A |vI/ ) = -f^) - I f) -If}... , (76) 



2 v y 3 
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where the solid line is the Green's function Gk(T2 — t\) and the vertex is 
Wide' (0(1") — 0(t — A)). An integration with respect to fictious time r is im- 
plied in each vertex and summation with respect to k in each line. 
The causal Green's function is given by the sum 

Gk = G fl f?(e k ) + G A 9(-^) 

where G R (G A ) are the retarded and advanced Green's function respectively 

Gk = -ie- %tk(T2 - Tl) 0{T 2 - n) 

and 

G A = ie- l€k(T2 - Tl) 6{Ti - t 2 ). 

In the diagrams of second and higher orders the main contribution in the 
sum over k comes from the region near the Fermi surface. For example, in the 
second order diagram we have 

(77) 

since the causal Green's function is a sum of advanced and retarded ones. The 
diagrams with only advanced or only retarded lines are equal to zero because of 
incompatible conditions set out by the time 9- functions. In the diagrams with 
both retarded and advanced lines because WW is large only if — < jr all 
energies are close to each other and therefore close to the Fermi surface. This 
allows us to substitute S(E F ,t) = S(t) for S(E,t) in Eq. @. We notice that 



d d k r(jf)i/(j0) fde 



(2-rr) d ■jfd/2—i / 2 7T 



dn. 



As the Green's functions do not depend on direction of k the integration over 
n = k/|k| involves only indexes of scattering matrix in each vertex. This in- 
tegration can be replaced by summation over scattering channels if we use an 
appropriate basis for the scattering matrix. 

The first diagram in ((76|) requires special attention because the contributions 
from all the region under the Fermi surface is significant and the energy depen- 
dence of S(E, t) can not be neglected. We notice, however, that this diagram is 
proportional to the average absorbed energy, because 

= -iA(0|i/oo - flo|0) - -iX(E) (78) 

The average absorbed energy was computed in |ll[ and is given by (||). Re- 
defining the Green's function and the vertex in the diagram expansion ( f76[) in 
a more convenient way we get (po|). 

C Integrating of continued g 



The integral in (|46|) with g given by (|47|,|48|) splits into two parts: the surface 
integral 

"(5) =-^[ tvid^gg-'d^gg-^drdx (79) 
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and the volume integral 

T(g) = -j^ J e^trid^gg-^gg^dxgg-^drdxdy. (80) 

Firstly, we will deal with the surface integral. Because of the identity 
tr(d l j,gg~ 1 d fl gg~ 1 ) = 4tr(d z gg~ 1 dzgg~ 1 ) the expression under the integral is 
zero in the region where r > or r < — a where d z g — 0. Combining 

trid^gh^g^d^gh^g- 1 ) = tr {d^gg^d^gg- 1 ) + 

+ tY{d^hh- 1 d^,hh- 1 ) + 2tr(g- 1 d^hh- 1 ) (81) 



with Eq. (|44|) we get 



1 '° 



57T 

f-0 



— a 



Sl(g) = -— dr dxtT{d x S- 1 (x)S(x)d x S- 1 (x)S(x)) 



~ -^J dT J dxtTiSdxS-idJ&f- 1 ^)). (82) 

Notice, that d x f = —idgf due to antianalyticity of /. 

In Eq. ( |S(i|) the interior of the sphere is separated into three parts by surfaces 
t = — a and r = 0. Everywhere except in the middle part where < t < 
—a the differential form ^{dgg" 1 A dgg~ x A dgg~ l ) is zero because in those 
regions g depends on two variables only. In the remaining part we have g = 
S^ 1 (x,y)f(z,y) and the differential form is exact 

trKdgg- 1 ) 3 } = - tr{(S- l dS) 3 } + tv{(df f~ x ) 3 } 

+ ZdtridSS- 1 Adff- 1 ) = Mtr{dSS- 1 Adff- 1 ). (83) 

Applying Stokes' theorem we get 

T(g) = ~ J dr J dxtriSdxS^d.fiz)/- 1 )- 

- JL f tYidSS^dfi-a + O-tx^y)/- 1 }- 

- -1- / tYidSS^dfi-O-ix^y)/- 1 } (84) 

where (3\ is a surface r = —a and @2 is a surface r = 0. The first term in ( |34| ) 
cancels the second one from (|2|) and thus we obtain (p9[). 
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